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Abstract: A procedure is reported for numerical analysis of false vacuum transition in a 
model with multiple scalar fields. It is a refined version of the approach by Konstandin and 
Huber. The alteration makes it possible to tackle a class of problems that was difficult or 
unsolvable with the original method, i.e. those with a distant or nonexistent true vacuum. 
An example with an unbounded-from-below direction is presented. 



This is a note on numerical computation of the decay rate of a false vacuum within 
a potential of multiple scalar fields. The aim is to describe a new strategy to obtain the 
bounce configuration, that is based on an earlier idea by Konstandin and Huber splitting 
the task into two stages [1]. The first is to search for an intermediate solution in one- 
dimensional spacetime using an improved potential. In the second, the field profile is 
stretched smoothly to that in four dimensions. 

One advantage of the original approach is that the outcome, by construction, must 
be either a valid solution of the Euclidean equation of motion or something clearly wrong. 
For instance, this property is not shared by the improved action method [2]. The authors 
of Ref. [3] use this method to find a fit and then check its quality by taking the kinetic to 
potential ratio within the action. This ratio must be —2 for a field configuration to leave 
the action stationary under an infinitesimal scale transformation [4, 5]. Being a necessary 
condition, however, this is not more than a consistency check and one cannot be sure that 
the numerical data is indeed a bounce until it is put into the equation. 

Its virtues notwithstanding, the original approach has drawbacks due to the use of an 
improved potential that render it suboptimal in some circumstances. Particularly prob- 
lematic is a case in which the false and the true vacua are separated far apart. This can 
arise for instance in the Minimal Supersymmetric Standard Model (MSSM) with a large 
trilinear term unless it involves a top squark [6, 7]. One can put an upper limit on such 
a term by demanding that the lifetime of the standard vacuum be longer than the age of 
the universe [8, 9]. The same criterion has been employed in a study of maximal flavour 
violation by a soft trilinear coupling [10]. For this last work, a modification has been made 
to the original approach, which is to be the subject of this article. 

Consider the potential depicted in Fig. 1, 
where it has been inverted to make it easy to 
imagine a classical motion. It is a function of 
scalar fields collectively denoted by <j>. The prob- 
lem is to estimate the rate of transition from the 
local minimum <f> + to the state with the low- 
est energy. Note that the notations <p± in this 
article follow Ref. [11] and are reversed with re- 
spect to Ref. [1]. 

Using a semiclassical approximation [11, 12], 
one finds that the decay rate of a false vacuum 
per unit volume is 

T/V = A exp(-SM), (1) 

where the determinantal factor A comes from 
the Gaussian functional integral around the stationary point eft, called the bounce, of the 
Euclidean action S. The formula for A can be found in Ref. [12]. It is very difficult to 
compute in practice and so can only be estimated to be around the characteristic mass scale 
of the problem. Normally, evaluation of A is exponentially less important than that of S[<f>]. 
Without loss of generality, the dominant bounce can be assumed to be 0(4)-invariant [5]. 
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Figure 1: An example potential V and the 
corresponding improved potentials U e and 
U e , displayed upside down. The false and 
the true vacua arc designated by 0+ and 
</>_, respectively. Between them, there are 
4> e and (ph. The former denotes any point 
with V(4> e ) = V (</>+) and the latter such a 
point on the bounce. 
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Therefore, one can let <j> be a function of the single coordinate p = \Jt 2 + x 2 , and write S 
in the form, 



oo 

2 / J„ „3 



S[</>(p)]=2ic* dpp 
J o 



(2) 



The ultimate goal is then to find the solution <ft of the Euler-Lagrange equation, 

W(0), (3) 
derived from the above Euclidean action, that meets the boundary conditions, 



d 2 cf> 3 dcf) 
dp 2 p dp 



dp 



(p = 0) = 0, cf)(p^ oo) = (f> + . 



(4) 



The original approach [1] is based on the idea that one can generalise the equation of 
motion to that in a spacetime dimensions, 



d 2 cf> a — 1 



+ 



dp 2 ' p dp 
which corresponds to the generalised action, 

/•oo 
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W(0), 



2 \dp 



+ v(4>) 



(5) 



(6) 



The procedure can be divided into two parts: to solve the equation in the undamped case 
with a = 1 and then to deform the solution to its final shape. The deformation part has 
no problem and shall be presented later when the refined version is described. 

The first part for the undamped case is summarised. If a = 1, the damping term 
in (5) disappears and so the motion conserves energy. Therefore, one can reproduce the 
initial condition in (4) in the limit where the starting and the ending points are degenerate. 
The degeneracy is approximated by flattening out part of the potential that is lower than 
the false vacuum. The solution is found by minimising the action (6) for a = 1. More 
specifically, the procedure is made up of the following steps. 



1. Find the minimum of Si with V in (6) replaced by the improved potential, 



v(<p) - 

Ue{9) = X + 



satisfying the boundary conditions, 



TO) - V{<t> + )] 2 | c2 



1/2 



0(O) = ^_, <f>(T) = </>+. 



(7) 



(8) 



The shape of U e is illustrated in Fig. 1 with e set to [V(<f)-) — V (<{>+)]/ 10. One can 
accelerate the convergence by making a further modification to U e to get 



U e (<f>) = U e {<j>) + AU € (<i>), 



0) 
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shown in Fig. 1, with a small perturbation, 



AW) = ~26 # + 3 e ,f"t +l l (10) 



• + l 3 !</>-- </>+| 2 
The initial choice of the parameter e and the time T may be 

emit = VtfJ) - V(<f>+), (11) 

r init = 20 1 0--<^+ i/Vw, (12) 

where Vj, is the height of the potential barrier. 

2. Find the point 4> e between 0_ and 4>+ such that V(4> e ) = V((j)+) within the configura- 
tion obtained in the previous step. Truncate the part before (fi e , i.e. the configuration 
on the plateau. 

3. Repeat minimisation with the new set of boundary conditions, 

0(0) = e , <(>(T) = </>+, (13) 

while iteratively sending e down to zero. Initially <p e is set to the point found above 
and then it is allowed to move freely on the submanifold with V{<p e ) = V(4>+). 

4. The last minimum for e = is the solution of the undamped equation, 

5? -ww, in, 

which is (5) for a = 1, and obeys the boundary conditions, 

^(0) = 0, <f>(T) = </,+. (15) 

There are a couple of points of concern. First, the improved potential (7), sketched in 
Fig. 1, has a plateau on which the solution of step 1 has to spend a long time. In fact this 
is what the improved potential is designed for so that the initial condition in (4) can be 
approximated. As already stated in Ref. [1], however, this costs generically many lattice 
points to store a long part of the path that is to be discarded eventually. This grows more 
and more problematic as 4>- moves away from fa, and hampers the computation. Second, 
step 3 might lead to a numerical instability. The intention of this step is to find a path 
that starts from fa and then rolls down the inverted potential. Although modified to have 
a plateau, the improved potential still has a part with U £ ((p) < U e ((f) e ) = £/ e (^>+) for a finite 
e. After starting from <p e , the path may prefer to climb up the inverted potential and spend 
time in this region to minimise the action. This can be avoided by using U e instead of U e . 
With a judicious choice of AU e , one can make 0+ the global minimum of U e so that the 
system tries to spend as much time as possible around (/>+ at the end of the path. Note 
that for this, e in (10) should be different from e in (7) in general. The initial condition 
in (15) is violated by U e (<j) e ) — V(0 e ), which eventually goes away as e — > 0. 
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The first point makes it difficult to analyse many interesting problems. For instance, 
vacuum transition triggered by the stau trilinear coupling in the MSSM is studied in Ref. [9]. 
The charge-breaking global minimum in the scalar potential is far from the local one due 
to the small tau Yukawa coupling. To overcome this problem, they first found the bounce 
of a temporary potential with two nearly degenerate minima using the above steps, and 
then made a continuation to the actual potential by iteration. In the course of deforming 
the potential, one should be careful not to introduce a singular behaviour. 

In what follows, a more streamlined procedure shall be presented. It does not cost extra 
lattice points to be truncated in the end. It does not need a deformation of the potential. 
It does not rely on the location of the true vacuum. Consequently, it can deal with a 
potential that has widely separated local and global minima or that is even unbounded 
from below. 

The main idea is to exploit the energy conservation in the undamped case, which was 
already mentioned in Ref. [1]. This feature enables one to replace the Neumann boundary 
condition of the original problem by a constraint on the potential. As pointed out above, 
fixing the potential at p = is not enough to have the desired solution be a minimum since 
a path can lower the action further by shooting up the inverted potential and staying in a 
region with V((f>) < V(<p+). The refinement here is to eliminate those paths by additional 
constraints. 

Including the preparation and the continuation stages as well, the new series of steps 
would be as follows. 

1. Find a point (f> e such that V((j) e ) = V ((/)+) that is on the other side of the barrier. 
For instance, one can walk over the barrier along a valley starting from <j) + until the 
level comes down back to V (</>+). 

2. Construct an initial configuration 4>{p) such that 



The simplest choice would be a step-like profile in which (j)(p) is constant except for 
one jump somewhere in the middle, possibly around p = 0.1T. 

Choosing a suitable T is important for efficient evaluation. Recall that (17a) can 
substitute for (15) only in the limit of energy conservation. Realising this limit 
would require infinite T, whereas a small T saves the computation time and storage. 
However, it is difficult to estimate an optimal T before doing any minimisation. 
Initially, one may try (12) with <p e instead of <j>— and Vb from the previous step, and 
then T can be adjusted later as explained below. 

Optionally to prepare a better initial profile, one could perform a single-field min- 
imisation of Si using as the field the position on the segment connecting and (fi e . 
One should repeat this process until T is adjusted appropriately. The criterion for 
this and the constraints on the potential are the same as in the next step. 

3. Find a minimum of Si obeying the boundary conditions, 



<t>(o) = <Ae, 4>{T) = 4>+. 



(16) 





(17a) 
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in combination with 



V{(j){p)) > V(</>+) for < p < T. (17b) 

The inequality (17b) should be enforced at every lattice point in order to prevent the 
unwanted 'upshooting'. 

One should repeat the minimisation with an increased T until energy is sufficiently 
conserved. It usually works to require that the average difference between the ki- 
netic and the potential energy densities be less than around 1% of the barrier height 
measured within the resulting profile. 

4. The minimum is a solution of the equation, 

with the damping term killed by a = 1, that complies with the boundary condi- 
tions (15). In comparison to (5), there is an additional offset parameter Ap that has 
been introduced to avoid a pathological behaviour near to p = [1] . 

5. Make a continuation to the damped case by gradually increasing a from 1 to 4 with 
Ap fixed around 0.75T, and then send Ap also to zero. One can choose to stop at 
a = 3 for tunnelling in a finite temperature system [13]. For each pair of a and Ap, 
one can linearise (18) using the series expansion of the right hand side, 

W(0) « W(0) + (<f>- 0) • VW(0), (19) 

and then iteratively solve it by matrix inversion. Boundary conditions are set by 
(15). 

6. The final solution is the bounce configuration. 

It should be reminded that steps 4 and 5 have been copied from Ref. [1]. Modifications 
have been made only to the steps for the undamped case. 

The proof of step 3 follows directly from that of step 1 in the original method [1]. First 
of all, consider a solution determined by the improved potential (7) and the boundary 
conditions (8). It is a minimum of the action since 0_ is the global minimum of U e . Next, 
take the limit of e — > 0. Since the path on the plateau makes no contribution to the 
action, the rest of the path becomes a minimum of Si subject to the constraints (17). This 
motion off the plateau satisfies the equation (14) and the boundary conditions (15) which 
are equivalent to (17a) by virtue of energy conservation. Note that the improved potential 
is invoked only for the proof but never appears in the practical numerical calculation. 

There is a subtlety regarding the bounce as a minimum of the action under the condi- 
tions (17). In the limit of T — > oo, there is a flat direction of minima that initially spend 
different amounts of time staying at (f) e before rolling down. Any of this family of minima 
qualifies in principle as a solution of the undamped equation. In practice, however, it is 
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desirable to find a path that does not have a long constant part to save computation re- 
sources. This is why the initial profile in step 2 was set to have a jump at a position close 
to p = 0. For a finite T, the flat direction is slightly lifted. The minimum search routine 
should be able to deal with this kind of continuously (quasi-)degenerate minima. Also, the 
termination condition should be tuned in such a way that the routine stops after reaching 
a point that is acceptable to step 4 even though it is not the exact minimum. This needs 
some experiments. 

The minimisation problem in step 3 is a nontrivial task termed nonlinear programming 
[14] which is a broad research area on its own. One possible strategy is the following (see 
e.g. [15]). An equality constraint is taken care of by employing a Lagrange multiplier. 
An inequality constraint is replaced with an equality by introducing a slack variable with 
a limited range. This limit is fulfilled by adding a barrier term to the function to be 
minimised. This term is a product of the barrier parameter fi and a barrier function of 
the slack variable. The solution of the original problem is obtained by solving a sequence 
of subproblems with fi successively decreasing towards zero. At each 'outer iteration', the 
subproblem with the fixed /i is solved by a damped Newton's method which involves 'inner 
iterations', and then this and the previous outcomes are used for making the initial guess 
for the next subproblem with a smaller [i. The method in Ref. [15] needs the 'outer loop' to 
be repeated several times. If the algorithm is too hard to implement, one can let a library 
do the job. For instance, the example exhibited below uses the Ipopt package [16]. 

It might be of interest to compare this and the original approaches for a problem 
that both can solve. Suppose that one introduces n lattice points each holding / unfixed 
field variables between + and 4> e (see Fig. 1). The constrained minimisation problem in 
step 3 roughly corresponds to applying Newton's method to a function of n(f + 2) — 1 
free variables, i.e. nf field values, n Lagrange multipliers, and n — 1 slack variables. This 
procedure is repeated, each time with a decreased barrier parameter. To use the original 
approach on the other hand, one needs additional N points from <p e to (ft- except <p e that 
has already been counted. The function to be minimised depends on (n + N)f variables in 
the first step and nf + 1 in the second including one Lagrange multiplier. Newton's method 
is repeated while iteratively driving e — > 0. Given only this information, it is difficult to 
tell which way is more efficient on general grounds. Nevertheless, it is obvious that a case 
with large N is better handled with the new strategy since each Newton iteration becomes 
very expensive with the original. This happens when 0_ is distant from 4> e . Another limit 
arises with large / that may favour the new approach over the original. In this case, the 
number of auxiliary degrees of freedom 2n — 1 in the former becomes much smaller than 
Nf in the latter. 

Now, the refined procedure is applied to a problem with two real fields as a demon- 
stration. The scalar potential is chosen to be 

V = ^ + ^ - + a\4>l - (20) 

This potential is unbounded from below. In this sense, it is qualitatively different from 
the examples in Ref. [1]. Clearly, one cannot impose the boundary conditions (8) since 
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Figure 2: Contours of the scalar potential (20). The local minimum is marked with <fi + . The 
white curve is the set of points 4> e at the same level as </>+. The clotted and the solid curves show 
the bounces in the undamped and the damped cases, respectively. 

there exists no (ft-. Therefore, the original method is not applicable unless the potential is 
altered. There is a minimum with V = at (f>i = 02 = that is local as long as A ^ 0. 
One can see the metastability easily by taking the field direction, (j>i = <f>2 = a, that cuts 
out the quartic term leaving 



Along this direction, a barrier starts from the origin and ends at the position where A a = 
mf + m\ after which the potential drops below zero and keeps falling down. 

The shape of the potential is illustrated in Fig. 2. The masses and couplings in (20) 
are set to be 



There is a force that drives the path off the straight directions with (pi = ±02, arising from 
the gradients of both the trilinear and the mass terms. A large difference between and 
m\ has been introduced to emphasise this dynamics. 

The white contour corresponds to the point e in Fig. 1. Finding any point on the 
curve is enough to complete step 1. Successful execution up to step 3 should lead to 
the undamped solution for a = 1 shown in Fig. 2. Notice that the starting point of the 
trajectory is on the white curve as enforced by (17a). Then one can proceed to take steps 
4 and 5 to deform this path to the final bounce for a = 4, Ap = 0. It is also partly visible 
in the same figure. 

The solutions are plotted as functions of the time in Fig. 3. As a increases, the starting 
point retreats from the false vacuum due to the damping term in (18). Notice that the 
path does not diverge indefinitely even with an unbounded- from-below potential. This 
is understandable from the fact that the original problem formulated in (3) and (4) has 
nothing to do with any information on the global minimum. For the same reason, the 
procedure described in this article does not need to explore deep into the potential for 
an unreachable bottom. One can put this solution back into the Euclidean action (2) to 





m\ = 3.0, m\ = 



.5, A = 1.5, 5 = 0.4. 



(22) 
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calculate the final answer, 



S[cj)} = 488. 



(23) 



In this problem, the potential has a symmetry with respect to the reflection fa —> —<p2- 
This pairs a bounce with the other that makes exactly the same contribution to the decay 
rate. Therefore, one must multiply the decay rate (1) by 2. 

It is straightforward to write the lat- 
ticised versions of the fields, the Euclidean 
action, the Euler-Lagrange equation, and 
the boundary conditions, and they can be 
found in Ref. [1] for instance. The same 
reference reports an analysis of discretisa- 
tion error. It is determined by the final 
iteration step for continuation to four di- 
mensions that is common to the present 
and the original methods. 
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